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Abstract— The power output and thermal efficiency of a finite-time, optimized, extra-terrestrial, 
solar-radiant Stirling heat engine have been studied. The thermodynamic model adopted is a regen¬ 
erative gas Stirling cycle coupled to a heat source and heat sink by radiant heat transfer. Both the 
heat source and sink are assumed to have infinite heat-capacity rates. Expressions are obtained for 
optimum power and efficiency at optimum power for a cycle based on higher and lower tempera¬ 
ture bounds. 


INTRODUCTION 

Since its conception in the early 1800s, the Stirling engine has periodically enchanted engineers and 
physicists because of the theoretical potential for approaching the Carnot efficiency. 1 - 2 Today, interest 
in the Stirling engine is again on the rise. 3 Among the reasons for this are great advances in materials 
technology, inherent environmental advantages of the Stirling engine and the fact that, as an externally- 
heated engine, it can be powered by a number of means. 4 The Stirling engine is either currently being 
used or is under consideration for use in a wide variety of applications, 5-8 including both extra-terrestrial 
and terrestrial solar installations. 9 - 10 

For lack of an alternative procedure, practical engineers have typically made use of the thermodyn¬ 
amic air-standard approach to model and study trends in applications of real gas Stirling engines. This 
analysis provides exceedingly generous estimates of the potential performance of such engines and is 
not relevant in the design activity. Thus, to provide a more reasonable estimate of the performance 
potential of a real regenerative solar-radiant Stirling engine, the innovative finite-time thermodynamics 11 
approach is employed herein on an endoreversible version of the Stirling cycle. An endoreversible cycle 
is one in which the external heat-transfer processes are considered the only irreversible processes within 
the cycle. 12 The preferred use of finite-time thermodynamics in this study is based on the inherent 
reality that a heat engine must produce work in a limited amount of time. In order to accomplish this 
task, the external heat-transfer processes of the cycle are modeled as occurring across finite-temperature 
differences and are thereby considered irreversible. The endoreversible cycle is then optimized with 
respect to power rather than efficiency (the usual focus of thermodynamic analysis). For clarity, this 
study is thus based on the recognition that real engines do not produce work over an infinite time but 
must produce power. Work produced over an infinite amount of time gives zero power. 

The relationships developed and the analysis obtained should provide insight into the operational 
and design criterion necessary for the attainment of optimum power in a simple solar-radiant Stirling 
engine for use in space applications. For this baseline study, only black-body surfaces are considered 
in the external heat-exchange process by radiation, The discussion is limited to a thermodynamic analy¬ 
sis of the engine itself and thus does not include a detailed description of the actual hardware involved 
in this energy-conversion process. 

THE ENDOREVERSIBLE EXTRA-TERRESTRIAL, SOLAR, RADIANT STIRLING ENGINE AND ITS FINITE¬ 
TIME POWER OPTIMIZATION 

The regenerative endoreversible Stirling cycle is depicted in Figs. 1 and 2. This cycle approximates 
the compression stroke of the real engine as an isothermal process 1-2, with irreversible heat rejection 
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Fig. 1. Functional schematic of the regenerative Stirling engine; note that the opposing pistons are 90 degrees 

out of phase. 



Fig. 2, Temperature-entropy diagram of the regenerative Stirling engine (ideal). 
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to an intermediate low-temperature sink. The heat addition to the working fluid from the regenerator 
is modeled as a reversible constant-volume process 2-3. The work-producing expansion stroke is mod¬ 
eled as an isothermal process 3-4, with irreversible heat addition from an intermediate high-temperature 
source. Finally, heat rejection to the regenerator is modeled using a reversible constant-volume process 
4-1. 

Within real Stirling engines, the external heat-transfer processes 1-2 and 3-4 must each occur in finite 
times. This in turn requires that these heat-transfer processes must proceed through finite-temperature 
differences and are therefore defined as being externally irreversible. During the first of these, the heat 
rejection process 1-2, energy flows from the working fluid which is maintained at a constant temperature 
r c , to the intermediate low temperature sink T L , thus flowing across the temperature difference (T c - T L ), 
as shown in Fig. 1. Similarly, during the heat addition process 3-4, heat is transferred from the inter¬ 
mediate high temperature sources T H to the working fluid which is maintained at a constant temperature 
T w , thus flowing across the temperature difference ( T H - 7 W ). 

For satellites in our solar system, the ultimate heat source is the sun (at temperature T s ) and the 
ultimate heat sink is outer space (at temperature T 0 ). If T H is the temperature at the external (radiative) 
surface of the black-body heater plate, then heat is transferred to this surface across the temperature 
difference (T s - T n ) and from this surface across (T H - T 0 ). Similarly, if T L is the temperature of the 
external surface of the cooler plate, heat transfer to the universe from this surface takes place across 
the temperature difference (T L - T Q ). 

If the regeneration is assumed to be ideal, the heat rejected to the regenerator by the hot working 
fluid during the process 4-1 is equivalent to the heat supplied to the cold working fluid during 2-3. 
This can be visualized by noting that the area under the process path on the temperature-entropy 
diagram of Fig. 2 for process 4-1, area c-l-4-d-c, is equal in magnitude to the area under process 
path 2-3, area a-2-3-b-a. The reasonableness of this validated by the fact that assumption is the 
efficiency of regenerators is continuing to improve, with at least one manufacturer as far back as 1972 
reporting regenerator effectiveness values of 95%, 13 and progressing to reported values of 98-99%. 14 
Thus, the relation for the overall ideal regenerative heat-transfer processes is 

2 G3 = 4G.- (1) 

Since the regenerator supplies the thermal energy transferred to the working fluid during 2-3, the 
only source of external heat addition to the cycle is provided by the intermediate source during process 
3-4. If the working fluid is considered to be an ideal gas, the magnitude of heat addition during 3-4 
can be found by substituting the definition for boundary work into the conservation of energy 
relation, yielding 


Gin = 3 G 4 = mRT v (\nr y ), 


( 2 ) 


where r v = compression ratio = V X IV 2 - Similarly, the external heat rejection 1-2 is 

Qo ui = - 1 G 2 = mRT c ( Inr v ). (3) 

For the thermodynamic cycle 1 -2-3-4-1, conservation of energy yields 

W„e, = Gnet = Gin " Gou, • (4) 

Substituting Eqs. (2) and (3) into Eq. (4) gives 

W na = mR(lnr v )(T w - T c ). (5) 

The average net power for a cycle is found by dividing the net cycle work by the time for one cycle 
t. This approach yields the following equation for the net cycle power: 


P = W„Jt = mR(lnr v )(T w - T c )!t. 


(6) 
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The total cycle time is the sum of the individual process times, i.e. 

* = *12 + *23 + *34 + *41 • (7) 

An order of magnitude analysis 4 shows that the times required to accomplish the isothermal com¬ 
pression and expansion processes are of the same order of magnitude as the times required to accomplish 
the constant volume processes; therefore all four times in Eq. (7) must be evaluated. 

Using Eqs. (2) and (3), the average rates of external heat transfer into and out of the engine are 
given by 


Gm = Cin /*34 = mRT w (lnr v )/t M (8) 

and 

Gout = Qou t /*i 2 = mRT c (lar v )/t n . (9) 

The rates of external heat transfer into and out of the engine may further be quantified from heat 
transfer theory as being proportional to the temperature between the respective intermediate thermal 
reservoirs and the constant temperature states of the working fluid during each process. One more set 
of expressions for these rates can be found from quantifying the black body exchange between these 
intermediate reservoirs and their respective ultimate thermal reservoirs. All of this implies that Eqs. (8) 
and (9) can also be expressed as 

Gm = UMT h - T w ) = A h a[F Hs n + (1 - F HS )T% - 7* H ] = A H o<F HS U ~ H) (10) 


and 


Qo Ut = U<Ac(T c - T l ) = A L F L o ain-r o ), (11) 

where, on the thermal side internal to the engine assembly, U represents an average overall heat-transfer 
coefficient due to conduction and convection between a radiative surface and the cycle working subst¬ 
ance and where the effective areas, A w and A c , represent heat transfer areas appropriate to the proper 
specification of thermal energy exchange between these radiative surfaces and the cycle working subst¬ 
ance during the respective warm and cold constant temperature processes. From the radiation side, A H 
and A l are the system areas receiving and emitting net radiant thermal energy, respectively. Also, u is 
the Stefan-Boltzmann constant and F HS and F LO are the view factors defined as the fractions of the 
energy leaving the surfaces represented by the first subscripts and reaching the thermal reservoirs rep¬ 
resented by the second subscripts. 

Substituting Eqs. (10) and (11) into Eqs. (8) and (9), respectively, yields expressions for the times 
associated with the external heat-transfer processes, viz., 

* 12 = mRr c (\nr v )/[U c A c (T c -T)] = mRT c (\nr y )/[A L F LO o{r L - T* 0 )] , (12) 

and 

*34 = m7?7’ w (lnr v )/[(/ w A w (7’ H - T w )] = m7?r w (lnr v )/[A H a(FHs7| - 7^)]. (13) 

Equation (6) yields expressions for the power which are based upon reservoir and isothermal working 
fluid temperatures. 

From Ref. 4, the times associated with the regenerative heat-transfer processes can be expressed as 

*23 = * 4 . = rnc v (T w - r c )/[C/ reg A reg ( 1 K)]. (14) 

Equations (12)—(14) may then be substituted in turn into Eq. (6), yielding expressions for the power 
which are based upon reservoir and isothermal working fluid temperatures, viz., 
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P = [TJ[UMTh -T„)(T w - T c )] + T C I[U C A C (T C - T L )(T W - TJ] 

+ 2c v /[/?(lnr v )t/ rcg A reg ]}-' (15a) 


and 


P = { TJ[A H o{F Hs n ~ r H )(T w - T c )] + T c /[A l F lo oOI - To)(T w - T c )] 

+ 2c v /[/?(lnr v )t/ reg A reg ]}-'. (15b) 


We note from Eq. (15) that power is a function of the variable temperatures F w , F H , T c , and 7’ L . The 
maximum power with respect to these four as yet undetermined temperatures can be obtained through 
the simultaneous solution of the following set of equations: 

dp/dT w = 0, 

(16a) 

o 

ii 

fO 

1- 

(16b) 

dp/dT H - 0, 

(16c) 

and 


dp/dT^ = 0, 

(16d) 

Application of Eqs. (16a) and (16b) to Eq. (15a) results 
working substance temperatures and the intermediate heat 
maximum power, namely. 

in functional relationships between the 
reservoir temperatures necessary for 

(F c )op, = CTt 5 

(17) 

and 


(F w ) op ,=cne , 

(18) 


where 


C = [( U^Th ) 05 + ( t/ATl) 0 5 ]/[( tVU 05 + ( UAc) 05 ) ■ (19) 

Proper substitution of these three results into Eq. (15b) enables the formulation for power to be recast 
into two new expressions, one in terms of the independent variables T H and T L only and the other in 
terms of the independent variables T c and 7" w only. The maximum power can then be obtained either 
by application of the Eqs. (16c) and (16d) to the first of these expressions or application of the 
Eqs. (16a) and (16b) to the second. To ensure the integrity of the work, both routes were employed 
yielding the same results for (T L ) opt and (T H ) opt . 

Assuming T 0 to be absolute zero, 


(TO = {7A h [F hs 71 - (7 h) op i] 2 /[A l F lo (F hs 7s + 7(7?,) opt )]) 1/4 (20) 


and 


(7-„) opt = {A H [F Hs n - (r„) opt ] 2 [F Hs n + 7(7' H ) 4 pt ] 7 /(7 7 A L F LO )} 1/36 • (21) 

Equation (21) is transcendental in nature and can be solved using a simple iterative numerical technique. 
The values given by Eqs. (17), (18), (20), and (21) had to be tested using a combination of techniques 
to ensure that they do in fact yield maximum (instead of minimum) values. Based on Eq. (15a) and 
these last two results, the expression for optimum power can be expressed as: 
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Pop, = {[[(t/A) 0 ' 5 + (tfvA)°' 5 ] 2 /[tfAtfwA( w, - (7 l)op,] 2 ] 

+ 2c v /[t/ reg A reg /?(lnr v )]}-'. (22) 

The work per cycle at maximum power or work optimum is determined by substituting the optimum 
working fluid temperature values from Eqs. (17) and (18) into Eq. (5), resulting in 

W opl = CmR(\nr v )[(T H )%] - (7^], (23) 

where C is defined in Eq. (19). 

Finally, the thermal efficiency of the cycle at the condition of optimum power is found by first 
substituting the optimum warm fluid temperature, Eq. (18), into Eq. (2) for the external heat transfer 
to the working fluid. This result and the optimum work given by Eq. (23) are then substituted into the 
basic definition for thermal efficiency, yielding 

T?„p, = W„ P ,/(fi„)opt = P(lnr v )[(7’ H )opi - (7’ I X 5 ,]/[(7’ H )2 p 5 t R ln(r v )] 

= 1-[(71)o P ,/(7’h)^, (24) 

Vupt = i - (t c /t w ) opi = i - [(r L /r H ) opt ]° 5 < l - r () /r H = Vth , (25) 

where 1 - (T Q /T H ) is the efficiency that could be obtained if the process were allowed to take place over 
infinite time in a cycle using an ideal regenerator and where (T H ) = T s (<tF h s ) 0 25 . Thus the efficiency at 
optimum power will always be a value that is substantially less than this upper bound. 

NUMERICAL EXAMPLE FOR A SOLAR SATELLITE 

Consider a satellite solar-radiant Stirling engine which receives radiant heat from the sun 
( T s = 5755 K) and emits radiant heat to space (T 0 = 0 K). The radius of the sun (R) is 6.95 x 10 8 m 
and the distance between the sun and the satellite (L) is 1.49 X 10 11 m. F HS = R 2 /(R 2 + L 2 ) = 2.176 x 10~ 5 
and F L o = 1- Thus the expression F HS a 1% contained in Eq. (10) has a power density of 1.353 kWm” 2 
which is known to be the solar constant. Values for A L and A H can be found through the use of Eqs. (12) 
and (13), respectively, once the other values contained in these expressions are specified. Note however, 
that these equations contain unknown temperature values. Yet, for initial computations A L and A H are 
not needed if time symmetry between the external heating process, 3-4, and the cooling processes, 1- 
2, is maintained. This symmetry condition requires that 


A h /A l = F h0 (Tl - r o m 5 /[(F Hs n - M 5 ] , (26) 

which can be used in turn to eliminate A H and A L from Eqs. (20) and (21). As noted in Ref. 14, time 
symmetry between these two processes also requires that (7 W A W = LA- To be comparable with typical 
existing Stirling engines, 4 values of C/wA w in the range 1-4 kW/K were considered for an engine with 
a regenerator U reg A reg value of 1000 kW/K, with a compression ratio (r v ) of 2 and with 9.3 x 10” 3 kg 
of air as the working fluid. A c v value of 0.7165 kJ/kg(K) and an R value of 0.287 kJ/kg(K) were used 
in the calculations. 

For the above installation, the values of (T H ) op „ (T w ) opl , (r c ) opt , and (TAp, were found to be the 
constant values 330.5, 293.8, 228.5 and 200.0 K, irrespective of Uj\„. Figure 3 contains a plot of P opl 
vs Uji w for this engine. The thermal efficiency of this power optimized engine is 0.222 for all L/„A W 
values, compared to a theoretical infinite-time upper-bound value of 1.0 for a Stirling engine operating 
between 5755 and 0 K and also using an ideal regenerator. 

A parametric study was conducted using this formulation to analyze the effect of F m on tj op , and 
Fop,. It was found even though increasing F HS by some means (such as by use of a parabolic reflector) 
can improve the power (Fig. 3), the value of i 7 opt remains unchanged. For a F HS of 8.0, for example, 
values for (TAp,, (T w ) op „ (T c ) opt and <T L ) op , were found to be 555.9,494.1, 384.3 and 336.3 K, respect¬ 
ively and the values for the power for different values of l^A were found to be approximately 1.682 
times those found for corresponding 1/*A W values the unenhanced F HS case (Fig. 3). However Tj opt is 
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Fig. 3. Optimum power vs heat conductance (U»A„) for various shape factors (F HS ). 


0.222 for this F HS for all values also. Thus one can see that an rj opl of 0.222 represents the upper 
limit for an engine with the above internal characteristics, performing at optimum power and using a 
simple black-body collector. 


CONCLUSIONS 

This work demonstrates the utility of a finite-time analysis approach in both the study and optimiz¬ 
ation of an extra-terrestrial, solar-radiant Stirling heat engine. The efficiency at optimum power obtained 
for the model studied is much more realistic than that provided by traditional thermodynamic techniques 
involving infinitely large cycle times. Also, while the latter approach results in the production of zero 
power, the finite-time approach, given herein, predicts very reasonable values. 
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NOMENCLATURE 


A C A w = Effective heat-transfer surface areas 
between the endoreversible Stirling 
engine and the intermediate low- 
and high-temperature reservoirs, 
respectively 

A h Al = Heat-transfer surface area of the 
heater and cooler plates, 
respectively 

A reg = Heat-transfer surface area of the 
Stirling engine regenerator 
C = Heat-engine parameter defined by 
Eq. (19) 

c v = Specific heat at constant volume 
(kJ/kg-K) 

Ehs<E lo = Shape factors between the 

subscripted heat-transfer surfaces 
(A h and A l ) and the subscripted 
thermal reservoirs (S—sun, 

O = outer space) 

L = Distance between the Stirling 
engine and the high-temperature 
reservoir 

m = Mass of working fluid in the engine 
P,P opt = Net cycle power output and 

optimum net cycle power output, 
respectively 

Qm,iQ* = Isothermal heat transfer into the 
engine from the intermediate high 
temperature source 

£?out. i Qi = Isothermal heat transfer from or to 
the engine to or from the 
intermediate low temperature sink, 
respectively 

2 Qi = Heat transfer from the regenerator 
to the working fluid 

4 j2, = Heat transfer captured by the 

regenerator from the working fluid 

Q m = Rate of heat transfer into the 

engine from the high-temperature 
source 

fiou. = Rate of heat transfer into the 
engine from the low-temperature 
sink 

R = Gas constant; radius of the sun 
r = Compression ratio, VI/V2 
s = Specific entropy 


T c = Temperature of the working fluid 
during the isothermal external heat- 
rejection process to the intermediate 
low-temperature sink 

T h ,T l - Temperatures of the intermediate 
heat source and intermediate heat 
sink, respectively 

T 0 = Temperature of outer space (the 
low-temperature heat sink) 

T s = Temperature of the sun (high- 
temperature heat source) 

T w = Temperature of the working fluid 
during the isothermal external heat- 
input process from the intermediate 
high-temperature source 
t = Total cycle time 
r, 2 = Time required to accomplish the 
isothermal compression 
r 23 = Time required for the heat-transfer 
process from the regenerator to the 
working fluid 

t M = Time required to accomplish the 
isothermal power stroke 
r 4 , =Time required for the heat-transfer 
process from the working fluid to 
the regenerator 

(/ c ,(/ w = Overall heat-transfer coefficient 
between the engine and the 
intermediate low-temperature sink 
or the intermediate high- 
temperature source, respectively 
U ng = Overall heat-transfer coefficient to 
and from the regenerator and the 
working substance during the 
constant-volume processes 
v = Specific volume 

V X ,V 2 = Volume at bottom dead center and 
top dead center, respectively 
WnM.Wop, = Net work output for one cycle and 
net work output per cycle at the 
optimum power condition, 
respectively 

rjop„rj, h = cycle thermal efficiency at the 
optimum power condition and 
infinite time, respectively 
cr = Stefan-Boltzmann constant 



